Polymer induced phase coexistence in systems of lamellar phases 
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The effect of adding nonadsorbing polymer to a lamellar phase of surfactant bilayers is studied 
theoretically. We find that the polymer produces coexistence between two lamellar phases of different 
layer spacings. The coexistence region is a closed loop, as in the experiments of Ficheux et al. [J. 
de Physique II 5 823 (1995)]. Within our model the coexistence is driven by depletion. 
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I. INTRODUCTION 

Many surfactants in solution form bilayers, i.e., two 
back-to-back layers of surfactants molecules with their 
alkane-like tails in the centre and only their hydrophilic 
headgroups in contact with the surrounding water |Q. 
These bilayers have been extensively studied both be- 
cause of their diverse behaviour J^] and because they 
serve as simplified model systems for cell membranes [|3| . 
In solution, bilayers form a stack of parallel layers: a 
lamellar phase. Here, we study theoretically the effect 
of adding non-adsorbing polymer to an unbound lamel- 
lar phase In particular, we attempt to understand 
some recent experimental results |^||| . A schematic of the 
system is shown in Fig. 1. An unbound lamellar phase 
is one in which the interactions between the bilayers are 
purely repulsive. The polymer-bilayer interaction is also 
purely repulsive, leading to a depletion in the densit 
of polymer segments near the surface of a bilayer J" 
This depletion has an associated free energy cost which 
is reduced if two bilayers approach each other, then the 
two depletion layers overlap [p|-jTT|]. Thus the polymer- 
induces an attraction between the bilayers. The attrac- 
tion leads to phase separation: two lamellar phases with 
different spacings between the bilayers coexist. 



II. MODEL AND THEORY 

The lamellar plus polymer phase is described starting 
from the free energy of a pure lamellar phase which is 
well known, see Refs. |j],|l2| 0- Then semidilute polymer 
solutions are discussed and the confining effect of the bi- 
layers on a semidilute polymer solution is estimated. The 
confinement free energy of the polymer solution between 
the bilayers is then added to the free energy due to the 
direct interaction between the bilayers to produce a fi- 
nal free energy. This free energy is used to calculate the 
phase behaviour of the bilayers plus polymer system. 



A. The interactions between the bilayers 

We treat the lamellar phase as a stack of uniformly- 
spaced rigid plates; thus we can define a separation d 
between consecutive plates in the stack. This is the spac- 
ing between the surfaces of consecutive plates; the smec- 
tic period is d + do, where do is the (assumed constant) 
thickness of a bilayer. Of course, the bilayers are not 
completely rigid but they bend on length scales that are 
much larger than d or the size of a polymer coil M . We 
neglect any interactions between nonadjacent plates. 

In the absence of polymer the free energy F of a stack 
of n layers, each of area A, is taken to be 



+ A h cxp(-d/\ h ), 



(1) 



riA d 

The first of the two terms of Eq. ([!]) is the contribu- 
tion of electrostatic interactions and is long ranged while 
the second term is due to the so-called hydration forces 
]l|, ^2|Jl4| . The hydration interaction is much shorter 
ranged than the electrostatic interaction, its range is of 
the order of lnm. We will not derive Eq. (Q) here, it is 
a quite standard expression for the interaction of a pair 
of charged bilayers separated by water. The electrostatic 
term is derived from the Poisson-Boltzmann equation for 
two charged plates separated by a gas of the counteri- 
ons [p]Jl3|; without added salt to screen the electrostatic 
interactions between the bilayers. Moderate screening 
has been found not to change the phase behaviour we 
predict. A e is a positive constant that is proportional 
to the surface charge density on the plates. The second 
term in Eq. (^) is essentially an empirical description 
of the strong repulsion between two bilayers less than 
about lnm apart. This repulsion is, to a good approxi- 
mation, exponential. Although it is strong, Ah 3> A e /d 
for d — lnm, this repulsion is much shorter ranged than 
the electrostatic interaction; is small ~ 0.2-0.35nm 
P4P^| . Equation ([j]) does not include the contribution 
to the interlayer interaction of the Helfrich or undula- 
tion forces P,|i~5|], because these are dominated by the 
electrostatic repulsions jl3j if, as here, the electrostatic 
repulsions are not screened. 
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B. The polymer— bilayer interactions 

As the polymer is nonadsorbing the bilayer-polymer 
interaction is assumed to be a simple excluded volume 
interaction, i.e., each bilayer excludes polymer segments 
from a planar volume. The polymer will be in the semidi- 
lute regime, we therefore use a simple scaling description 
of the structure and thermodynamics of the polymer so- 
lution The size of an isolated polymer coil of N 
segments each of length a, is measured by its radius of 
gyration Rq which is given by Rg ~ aiV 3 / 5 , in a good sol- 
vent The concentration of the polymer is specified 
by its volume fraction fa. The subscript b is to remind 
us that it is the volume fraction of the polymer far from 
the surface of a bilayer. Near the surface of a bilayer the 
polymer's density drops below fa. The average volume 
fraction of polymer in the surfactant plus polymer sys- 
tem is not equal to fa] it will be less. The correlations 
(fluctuations) in an isolated polymer coil extend across 
the entire coil; the correlation length £ is thus £ ~ Rg- In 
a bulk solution of polymer, when the polymer coils begin 
to overlap, interactions between different chains reduce 
the correlation length £ ]l6| . The overlap volume fraction 
4>b * s 01 = N~ 4 / 5 . This volume fraction is the boundary 
between the dilute and semidilute regimes; below fa the 
interactions between chains can be treated as a pertur- 
bation but above fa the interaction between a pair of 
chains is 3> ksT, where ks is Boltzmann's constant and 
T is the temperature. So, £ = Rg in the dilute regime, 
and 



-3/4 



in the semidilute regime. (2) 



solution is no longer uniform, the polymer density at a 
point will depend on the distance from the bilayers above 
and below the point. However, at distances ^> £ from any 
bilayer, the polymer density will reach a constant value. 

We model the bilayer-polymer interaction as a hard 
wall-polymer interaction [7|; there is no adsorption. The 
correlation length £ is a measure of the distance over 
which correlations decay. Therefore, if we introduce a 
single bilayer, modelled as a hard sheet, into a polymer 
solution it perturbs the solution above and below its sur- 
face out to a distance of order £. This is a much longer 
ranged perturbation than in a simple liquid where the 
density of the fluid returns to its bulk value in a distance 
of a few times a. Within £ of the surface the conforma- 
tional freedom of the polymer chains is reduced by the 
presence of the surface, reducing their entropy and hence 
their density near the surfaces of the bilayers 0-^] . The 
introduction of the bilayer has changed the free energy 
of the polymer solution by introducing two polymer-wall 
'interfaces' (they are not true interfaces as they do not 
separate two coexisting bulk phases). This layer of re- 
duced density is often called a depletion layer [[n]. We 
require the difference 7, between the free energy of unit 
area of the layer which is depleted in polymer and the 
free energy of unit area of a layer of the same width of 
the bulk fluid. It is, essentially, the free energy of a hard 
wall-polymer interface. The picture is of a layer of thick- 
ness £ in which the density drops from its bulk value fa 
to a much lower value. The excess free energy of the de- 
pletion layer is then approximately equal to the amount 
of work done in squeezing out the polymer from a slab of 
unit area and width £ JT^]. Thus, 



We require the form of the osmotic pressure of a polymer 
solution, it is 



(6) 



no 3 



in the semidilute regime. (3) 



The Helmholtz free energy of a polymer solution F p is 
also required. We can obtain the density dependent part 
of the free energy, which is all that is required for de- 
termining phase coexistence, from the thermodynamic 
relation II = —dF p /dV, where V is the volume of the so- 
lution. The volume derivative may be transformed into a 
density derivative. Then the integration of II as a func- 
tion of density leads to 



k B TV 



const. 



(4) 



Finally, the polymer's chemical potential fj, p is obtained 
from n p = F p a 3 /(Vfa) + Ua 3 /(f> 



~ const. + 4>l ,i . 

kb-1 



(5) 



All of the above is for a bulk polymer solution, i.e., a 
solution just of polymer. In the presence of bilayers the 



This is the surface tension of one wall-polymer solution 
interface. In a system of bilayers and polymer we can 
use this expression provided both II and £ correspond 
to those values in a bulk solution at the density of the 
polymer solution far, 3> £, from the surface of any bi- 
layer. This can perhaps be made more clear if we con- 
sider that the surfactant plus polymer system is in equi- 
librium with a pure polymer solution across a semiper- 
meable membrane permeable to the polymer but not to 
the surfactant. Then the volume fraction of the pure 
polymer solution in equilibrium with our bilayers plus 
polymer solution, is the volume fraction at which the II 
and £ of Eq. (g) should be evaluated. 

Each bilayer creates two polymer depletion layers. If 
the bilayers are sufficiently far apart, d ^> £, then these 
layers are independent and we have 2n depletion layers, 
each of which contributes an amount to the free energy 
given by Eq. (|^). However, as the bilayer separation d is 
reduced to ~ £ the depletion layers of adjacent bilayers 
overlap. When the layers overlap, the volume occupied 
by these layers is less and therefore the total free energy 
associated with these layers decreases. It is this reduction 
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in free energy which, we suggest, is responsible for creat- 
ing coexistence between two bilayer phases with different 
bilayer spacings. 

A difficulty now arises: the analysis which led to our 
estimate of 7 in Eq. (^J) was purely qualitative, whereas 
we now require a definite expression for how the free en- 
ergy of two depletion layers changes as they coalesce. The 
definite expression is required to calculate a theoretical 
phase diagram in order to compare with experimental 
phase behaviour. The free energy of a pair of depletion 
layers 72(d), must be 27 if there are far apart, d ^> £. 
Within the scaling approach there is only one relevant 
length scale: £. Thus, 72 is of the form 



72(d) = 2ne/(d/0; 



(7) 



that is, we are assuming that the density dependence is 
only through £. The functional form of /(d/£) is the 
same at all densities. When d <C £ the polymer density 
between the two bilayers will be much less than in the 
bulk at the same chemical potential. But if the polymer 
density between the plates is negligible then the pressure 
exerted by the polymer inside the plates is negligible and 
the pressure due to the polymer solution is the negative 
of the pressure of the bulk polymer solution, indepen- 
dent of the separation d. If the polymer contribution to 
the pressure is a constant then its contribution to the 
free energy must be lid, which forces /(d/£) = d/(2£), 
for d <C £. In the other limit, d 3> £ 72 = 27 and so 
/(d/0 = 1- A function which satisfies both these re- 
quirements is f(z) = tanh(z/2). This is merely an ad 
hoc function which is consistent with the little we defi- 
nitely know about the free energy of a polymer solution 
in a slit [pp7[]. Brooks and Cates |l^] have taken a dif- 
ferent approach to the same problem. They obtain a free 
energy expression valid for all values of d by forcing the 
different scaling theory free energies in the different 
scaling regimes to match. Polymer induced coexistence 
occurs for d/£ ~ 1 and so it is for these values of d/£ 
that f(z) is required. However, scaling theory is unable 
to predict the form of the free energy when the parame- 
ter d/£ is of order unity. So, where we find coexistence, 
there is no reason to expect the free energy expression of 
Brooks and Cates to be any more accurate than the much 
simpler interpolation formula given here. They found co- 
existing lamellar phases, as we will do, but they did not 
find reentrance. In order to find reentrance it is necessary 
to explicitly include the contribution of hydration forces, 
which Brooks and Cates did not do. 

For the bilayer plus polymer solution system we use 
the variables d and /i p , as these are most convenient. 
A thermodynamic potential T/(nA) for a pair of adja- 
cent interacting bilayers is then obtained by adding the 
Helmholtz free energy due to the repulsions between the 
bilayers F/(nA), to the grand potential of the polymer 
between the bilayers, where the grand potential includes 
an interfacial term, 72. We treat the slit between a pair 
of bilayers as a thermodynamic phase. The grand poten- 
tial of the polymer in the slit is then, by definition p9[, 



equal to — HAd + J2A, as the volume of the phase is Ad. 
As we are fixing the chemical potential not the density 
of the polymer the bulk part of the contribution of the 
polymer is not its Helmholtz free energy but — T1V. So, 
from Eqs. (0) and (|) 



r 

nA 



F 
nA 



-nd + 2ne/(d/o 



(8) 



The right hand side of Eq. (g) is T per slit. The three 
terms are the direct repulsion between the two bilayers, 
the bulk part of the grand potential of the 'phase' be- 
tween the slits and the interfacialpart of the grand po- 
tential of this 'phase'. Equation (||) is correct in the two 
limits of pure bilayers and pure polymer, and the third 
term is a crude approximation to the additional cost of 
the polymer-bilayer depletion layer. In the derivation of 
Eq. (0) we have implicitly treated the bilayers plus poly- 
mer mixture not as a mixture but as a polymer solution 
between a stack of hard walls Because of this assump- 
tion n is unaltered by phase separation and is the same 
in the coexisting phases. It is not the pressure of the 
bilayer plus polymer system. 

Coexistence can be found from the common tangent 
condition. This leads to the requirement that at coexis- 
tence the two phases have equal values of d[T/(nA)]/dd 
and T/(nA) - d(d[T/(nA)]/dd). These two equalities de- 
termine the values of d in the two coexisting phases. 



III. COMPARISON WITH EXPERIMENT 

Ficheux et al. |(| studied a lamellar phase of the surfac- 
tant sodium dodecylsulphate (SDS) in an aqueous solu- 
tion of the polymer poly(ethyleneglycol) (PEG). We set 
the parameters of T (||) to be reasonable for their system. 
There are 5 parameters: A e , Ah, Xh, a and Rg- For the 
purposes of calculation we express them in reduced units. 
These units are defined using the length of a polymer 
segment a as our length scale and as our energy scale 
we use ksT. Then A e , Ah, Xh and Rg are expressed in 
units of ksT/a, ksT/a 2 , a and a, respectively. We take 
ksT = 4 x 10~ 21 J and a — 0.3nm; the first is approxi- 
mately its value at room temperature and the second is 
reasonable for a flexible polymer. We set A e — 0.5fceT/a, 
Ah = 50fcsT/a 2 , Xh — a and Rg = 10a. These values 
correspond to A e = 6.67 x 10 _12 Jm _1 , A h = 2.22Jm~ 2 , 
A = 0.3nm and Rg = 3nm. Ficheux et al. report an 
Rg = 2.9nm ||. A hydration decay length of 0.3nm is 
approximately what has been found in a number of ex- 
periments |l|.|l2 14 1 . Our value of Ah is comparable to 
the estimate given by Roux and Safinya of 2Jm~ 2 |]J 
and our A e is comparable to the 4.49 x 10 _12 Jm~ 1 ob- 
tained from (the dominant term in) the solution of the 
Poisson-Boltzmann equation fl3f| . 

Using this set of parameters, we have calculated the 
phase diagram of a stack of bilayers in the presence of a 
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polymer solution; the II— d plane of the diagram is shown 
in Fig. 2. There is a closed region of phase coexistence, 
as in the experiment of Ficheux et al; it is bounded at 
high and low polymer pressures by two critical points. 
Note, however, that Ficheux et al do not rule out the 
possibility of the polymer adsorbing onto the bilayers. If 
too little or too much polymer is added no phase coex- 
istence is found. The presence of a closed loop of phase 
separation is distinctly unusual. We will show why the- 
ory predicts it in the hope that our theory provides a 
reasonable description of the experiment of Ficheux et 
al. It is worth noticing that the closed loop is entirely 
within the semidilute regime. The dashed curve in Fig. 
2 shows 2£ which only reaches its 2Rq, indicating the 
crossover to the dilute regime, at values of II below those 
of the phase coexistence loop. 

In the absence of polymer, the free energy (Q), and 
the pressure, of a stack of bilayers is a monotonically de- 
creasing function of d. The plates repel each other at 
all separations due to the repulsive electrostatic and hy- 
dration interactions. In the presence of polymer, there is 
an additional term in the thermodynamic potential (^), 
from the depletion layers. It is the overlap of these lay- 
ers that is driving phase separation. For a stack with 
d a few times £, the depletion layers contribute nearly 
27 per bilayer. However, if this phase separates into two 
phases: one with a large d and one with a d < 2£, then 
the total free energy of the polymer decreases because of 
the reduced total volume of the solution within £ of a 
bilayer surface. Thus the total contribution of the deple- 
tion layers to the free energy is reduced, at a cost in free 
energy due to the bilayer-bilayer interactions. As the 
parts of the free energy from bilayer-bilayer interactions 
increase more rapidly than linearly, the phase separation 
necessarily increases this part of the free energy. 

The total free energy available from overlapping the 
polymer depletion layers is of the order of Thus, 
phase separation will only occur if this is sufficient to 
overcome the repulsion of the bilayers. Of course, this 
repulsion decreases as d increases; so if £ is large the 
polymer depletion layers overlap when the bilayers are 
far apart and repel each other only weakly. Conversely, 
if £ is small the repulsions are strong. The polymer so- 
lution can only overcome the bilayer repulsions and pro- 
duce two bilayer phases if both IT and £ are sufficiently 
large. In the semidilute regime £ decreases as IT increases 
and so the two requirements are, to an extent, antago- 
nistic. They are only satisfied for a range of IT. 

The pressures due to the electrostatic and hydration 
interactions are A e /d 2 and (Ah/\h) exp(— d/Xh), respec- 
tively. In order for the polymer to control the phase 
behaviour the osmotic pressure of the polymer must be 
greater than both of these, at the separations at which 
the polymer depletion layers begin to overlap. So, we 
have two inequalities which must be satisfied in order for 
there to be phase separation, 

n > A e /e. (9) 



and 

n>^ e -e/^. (10) 

Aft 

Both of these inequalities only give the correct order of 
magnitudes. The correlation length £ varies with IT as 
£ = a(IT/T) -1 / 3 . Using this relation we see that the 
right hand side of Eq. (0) varies as IT 2 / 3 , more slowly 
than the left hand side. Thus Eq. (^) provides a lower 
bound to the polymer osmotic pressure required to induce 
phase coexistence. For values of IT above this bound the 
overlap of the interfaces releases enough free energy to 
overcome the electrostatic repulsions. Now, consider Eq. 
(|Io|). If we express £ in Eq. ( |i"o| ) in terms of IT we see 
that the right hand side of Eq. (|l0|) increases with IT 
as exp(— IT -1 / 3 ). This increase is more rapid than linear 
and so Eq. ( [To| ) provides an upper bound to the polymer 
osmotic pressure required to induce phase coexistence. 
The presence of these lower and upper bounds — due to 
the long and short ranged parts of the bilayer repulsions, 
respectively — leads straightforwardly to a closed loop 
coexistence region. There is slight subtlety in that the 
loop is d ~ 15 which is 5A^. The hydration repulsion 
is able to have an effect at such long range because the 
electrostatic repulsions and polymer induced attractions 
are delicately balanced in the coexistence region. Anal- 
ogous behaviour has been predicted by Cast et al. || 
when polymer is added to a suspension of charged spher- 
ical particles. There, as here, the decreasing range of 
the polymer induced attraction can cause reentrant mis- 
cibility. When the spheres repel each other with a long 
ranged electrostatic repulsion then the range of the poly- 
mer induced attraction between the spheres must have a 
longer ranged than this repulsion 

As polymer is added the dielectric constant of the so- 
lution between the bilayers decreases [^0| ; decreasing the 
repulsive electrostatic interactions. Our neglect of this 
effect may cause us to underestimate the size of the co- 
existence loop but incorporating a polymer concentration 
dependent A e will not change the qualitative phase be- 
haviour. 



IV. CONCLUSION 

Our theoretical model yields a closed-loop phase- 
coexistence region in the phase diagram. The model is a 
simple, additive, combination of a standard model for the 
interaction between charged bilayers and a crude scaling 
form for the free energy of a polymer solution in a slit. 
It provides a description of the behaviour of a mixture 
of surfactant bilayers and semidilute nonadsorbing poly- 
mer on a length scale of l-10nm; the length scale on 
which the entropy of the polymer solution and the repul- 
sive bilayer-bilayer interactions are competing. Thus, we 
provide an explanation of the phase behaviour in terms of 
behaviour at mesoscopic length scales. The model shows 
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qualitatively the same behaviour as found in experiment 
U, without any fitting. This encourages us to believe 
that our theory describes correctly the mechanism be- 
hind the experimental observations. 
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FIG. 1. Schematic of a pair of bilayers with a polymer solution between them. The thick curves denote the surfaces of the 
bilayers and the polymer is represented by the entangled thinner curves, d is the spacing between the bilayers and £ is the 
polymer correlation length. Note the low density of polymer near the surfaces of the bilayers: the depletion layer. 
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FIG. 2. The phase diagram of a bilayer plus polymer solution system, as calculated from the thermodynamic potential 
(@), with the parameter values given in the text. II is the osmotic pressure of the polymer and d is the spacing between 
the surfactant bilayers. The solid curve separates the one and two-phase regions and the dashed curve is 2£ as a function of 
the osmotic pressure of the polymer II. As II decreases £ increases until it reaches the radius of gyration, equal to 20d here, 
whereupon it remains constant. The thin horizontal lines are tie lines. 
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